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Experiments on spherical particles in a 3D Couette cell vibrated from below and sheared from 
above show a hysteretic freezing/melting transition. Under sufficient vibration a crystallized state 
is observed, which can be melted by sufficient shear. The critical line for this transition coincides 
with equal kinetic energies for vibration and shear. The force distribution is double-peaked in the 
crystalline state and single-peaked with an approximately exponential tail in the disordered state. 
A linear relation between pressure and volume (dP/dV > 0) exists for a continuum of partially 
and/or intermittently melted states over a range of parameters. 

PACS numbers: 45.70.-n, 64.70.-p, 64.60.Cn, 64.60.My 



The relative stability and selection of crystallized and 
disordered states is fundamental to the study of con- 
densed matter systems. For systems out of equilib- 
rium, the thermodynamic picture under which temper- 
ature melts crystalline order is not necessarily valid. 
Fluctuations provided by external driving, seemingly 
temperature-like, have been observed to either order or 
disorder a system. In particular, the mechanisms by 
which colloidal suspensions shear-melt and shear-order 
0, 01 , remain the subject of debate. Granular mate- 
rials, which are athermal and have strongly dissipative 
interactions, provide a complementary means to investi- 
gate the transition between disorder and crystallization 
in nonequilibrium systems. 

We examine a granular system in which there is a novel 
phase transition associated with competition between or- 
dering via one type of energy input, and disordering via 
another. We perform experiments in a variation on the 
Couette cell, a convenient setup for the study of sheared 
granular materials 0, S • An annular region containing 
monodisperse spheres is vibrated from below and sheared 
from above, with the mean pressure and volume set from 
below by a spring. In such a system, the shear and vibra- 
tion provide competing effects: the system is disordered 
or crystallized depending on their relative strengths, and 
the boundary between these two states occurs at equal 
kinetic energy input from these two driving mechanisms. 
The freezing transition is hysteretic, similar to a "freez- 
ing by heating" |^ transition. In addition, we observe 
metastable states over a range of packing fractions and 
forces. These states have the remarkable property that 
increased internal forces occur for less dense packings, in 
contrast to conventional mechanical or thermodynamic 
compressibility. While previous studies have investigated 
vibrated or continuously sheared granular materi- 
als separately, the combination of the two leads to novel 
effects. 

Fig. n shows a schematic of the apparatus. Monodis- 
perse polypropylene spheres fill the region between two 
concentric, stationary side walls, with a rotating upper 
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FIG. 1: Schematic cross-section of Gouette cell experiment 
(not to scale). Sample images at outer wall, showing crystal- 
lized (a) and disordered (b) phases, with glued black particles 
not visible on upper shearing plate. 

plate attached to a motor and a piston-like bottom plate 
attached to an electromagnetic shaker with a static in- 
ternal spring constant k = 341 N/m. The outer wall 
is Plexiglas to allow visualization and the bottom plate 
contains a force sensor flush with its smooth surface. The 
shear plate surface is a disordered layer of particles. The 
cell contains 71200 (0.43 kg) polypropylene spheres of 
diameter d = 2.337 mm and 5% polydispersity. 

The available parameter space is large: shearing rate Q 
for the upper ring, vibration amplitude A and frequency 
/, peak acceleration F = A(27r/)^/^ for sinusoidal vibra- 
tions, packing fraction 7, height and compressive load 
L. We set Q and F for fixed / (60 Hz) and number of 
particles; 7 varies freely. When the system is at its most 
compact, L is at its lowest; the cell is positioned so that 
this state has Lmin ^ (zero load) for F = 0. We explore 
shearing rates ofQ = 0.01 to 1.67 Hz (4.2 to 700 d/s) and 
peak accelerations F = to 7. Even for F = 7, the vi- 
bration amplitude A is only 0.2d^ small compared to the 
diameter of the particles. Relative humidity was 31% to 
35% for all runs; higher humidities shift the transitions 
observed without affecting the qualitative behavior. We 
measure forces on a circular plate (diameter 5 Ad) cen- 
tered on the bottom of the cell. Force measurements are 
done capacitively, to allow for high speed kHz) and 
high sensitivity (~ 10~^ N). The bottom plate moves in 
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FIG. 2: Hysteretic transition from disordered (right) to crys- 
tallized (left) state for F = 2.0. Dashed line is compaction 
achieved for = Hz, for comparison. The minimum 7 does 
not agree with jhcp = 0.74 due to the curved geometry and 
error in the measurement of ball and cell dimensions. Vertical 
bar represents standard deviation of height measurements. 



response to both the shaker and the particles. A laser 
position sensor measures the bottom plate position and 
a piezoelectric accelerometer measures F. Error in the 
bottom plate position is ±0.03 mm, and in H is 0.2 mm. 

The system crystallizes when sufficiently vibrated, 
with sustained vibration at F > 6 resulting in a hexago- 
nally close packed (HOP) state with a few (< 10) grain 
boundaries and point defects (holes) in the outer visible 
layer. We refer to this state as crystallized Exami- 
nation of the top surface suggests that the order persists 
throughout the whole cell. Less vigorous vibration re- 
sults in other HCP orientations and packings in some 
regions, as well as disordered regions. As seen in the up- 
per photo in Fig. ^ the crystallized state is sufficiently 
compact that only a few particles are in contact with the 
shearing surface. In fact, for higher loads, crystalliza- 
tion is suppressed since the Lmin ^ state is no longer 
accessible. 

In contrast, shearing with the upper plate has a disor- 
dering effect. The degree of disorder depends on the pa- 
rameter regime, and the effect is most pronounced near 
the top of the cell. Running the system at = 1.61 
Hz and moderate F = 2.0 prepares it in a fully disor- 
dered state. At the beginning of each run, we prepare 
such a state before setting the desired (l^,r). Follow- 
ing such a change, we observe transients of up to five 
hours. Changes in H are not necessarily monotonic and 
some states never settle to a characteristic mean an 
important effect discussed below. 

Fig. [2I shows the hysteretic behavior of transitions be- 
tween the disordered and crystallized states. As we de- 
crease Q in steps from the prepared disordered state, the 
decrease in H is logarithmic in Q . A first-order phase 
transition to the crystallized state occurs at ftc = 0-03 
Hz. Further decreases in Q result in only slight additional 
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FIG. 3: (a, b, d, e) Spatial autocorrelation functions (via 
Wiener-Khinchin method) of outer wall images for F = 2.0 
and decreasing Q starting from disordered state, (c) The av- 
erage magnitude of the four off-axis nearest neighbor peaks 
(x) of the autocorrelation functions. 

compaction. When we increase Q from this crystallized 
state, the transition to disorder is higher, at Qh = 0-08 
Hz. The shearing causes a transition to less compact 
state. The total change in H corresponds to 0.8(i, or 5% 
of the cell height. While this is a small change in 7, pre- 
vious 2D studies pjj indicate that a 4% decrease in 7 is 
sufficient to cross from a state with a strong force chain 
network to a fluid-like state without force chains. 

For F = 2.0, we measured the degree of order as a 
function of Q via 2D spatial autocorrelation functions of 
images obtained at the outer wall (see Fig.|2|). The mag- 
nitudes of the off-axis peaks, marked x in (b), change 
sharply at There is some order present even for 

Q > due to clusters, particularly in the bottom half 
of cell. Correlations in the direction parallel to the shear 
remain even as the transverse correlation disappears dur- 
ing the shear melting for Q > Qc- The remaining longi- 
tudinal correlations correspond to linear "trains" of par- 
ticles, and in separate experiments under larger confining 
pressures have been observed to be the dominant feature 
[3. The importance of clusters in the stability of granu- 
lar systems has been much investigated , and further 
work is necessary in this system. 

In addition to the visible changes in 7 and ordering, 
there are changes in the interparticle forces as well. Prob- 
ability distribution functions (PDFs) of the forces mea- 
sured at the bottom of the cell are markedly different 
for the two states. As shown in Fig. for Q < Qc 
the PDFs are double-peaked, corresponding to the two 
peaks in the PDF of a sinusoid. Thus, we infer that the 
crystallized particles respond to the motion of the bot- 
tom plate as a solid body. Time traces of force data 
with Q > Qc show amplitude disorder in the force re- 
sponse, but a well-defined frequency, /. The correspond- 
ing PDF is single-peaked and the tail of the distribution 
decays approximately exponentially, a shape characteris- 
tic of granular materials J-3]. The PDFs for increasing 
and decreasing Q show the same transition, but at Qh 
and rtc^ respectively. 

As measured from the kurtosis, the PDF is increasingly 
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FIG. 4: Force PDFs (scaled by mean) for representative 
decreasing-Q (<) points in Fig. |5] Solid lines are disordered 
and dashed are crystallized. Inset: Kurtosis of the PDFs, 
showing divergence at Qc Dashed line is a guide to the eye. 

wide 8iS Q ^ Qc (see Fig. 0J. This is reminiscent of the 
large fluctuations near a jamming transition in sheared 
colloids the broadening of the force PDF for simula- 
tions of sheared Lennard- Jones particles JjOj , and the 
divergence of the viscosity in classical glass transitions. 
While a crystallized state is fundamentally different from 
a jammed one due to the presence/absence of order, they 
share a lack of freedom to rearrange. 

As is common with sheared granular materials, the sys- 
tem exhibits an approximately exponential velocity pro- 
flle localized to the shearing wheel, with a decay length 
of a few ball diameters. This effect is most extreme in 
the crystallized state, where the system exhibits stick- 
slip behavior for the few balls in contact with the shearing 
wheel, similar to that observed in granular friction exper- 
iments 17J. Even in the disordered state, the vibrations 
do not create a more liquid-like state characterized by an 
approximately linear velocity proflle. 

The autocorrelation functions in Fig. |S1 show that dis- 
order occurs throughout the cell for Q > This full- 
depth disordering stands in contrast to the disordering 
localized to the shear surface in molecular dynamics sim- 
ulations of a 2D Couette cell with constant-P boundary 
conditions and gravity, but without vibration 18]. Since 
these simulations showed no disordering below the shear 
band, it appears that vibrations actually augment the 
disordering process. The effect we observe could be due 
to either the fluctuations induced by the vibration itself, 
or the periodically increased contact with the shearing 
surface. An initially crystallized state sheared without 
vibrations fails to disorder the lower portion of the cell, 
indicating that melting is suppressed. The vibrations 
therefore act subtly: promoting crystallization at low 
but also allowing low-shear areas of the cell to disorder 
at high Q. 

The crystallization/disordering transition can be un- 
derstood in terms of the relevant energy scales. Ignoring 



FIG. 5: Phase diagram in Q and F. Observations of crystal- 
lized (O) and disordered state (□). Overlapping circles and 
squares are hysteresis from FigEl Dotted line is En = Er. 

the effects of a compressive pressure {Lmin is small at 
l^c), shearing and vibration are the only two compet- 
ing effects. For each, there is an associated kinetic en- 
ergy input for a particle of mass m: Eq = ^m{RQ)'^ 
and Er = |m(27r/A)^. (These two modes are also 
anisotropic: the shear acts perpendicular to gravity and 
the vibration parallel.) For E^ > Eq the system is crys- 
tallized and for Er < Eq the system is disordered. The 
agreement of this line with experimental observations is 
shown in Fig. Hysteresis occurs because the pressure 
at the line Er = Eq depends on the height of the cell. 

The pressure P{t) on the force sensor shows signifl- 
cant fluctuations, but the average P(V), calculated by 
binning the data by volume V, is a smooth function (see 
Fig. Except at the highest ft and F, the data for 
non-crystallized states is consistent with a single value 
for dP/dV^ including during transients or in runs with 
higher Lrnin or higher humidity. At the highest values 
of 1], Fig.int shows a departure from the linear relation- 
ship (smaller dP/dV), coinciding with shorter tails in the 
high ft force PDFs. Signiflcantly, dP/dV is positive, cor- 
responding to a negative compressibility. 

In several experimental runs, the system intermittently 
exhibited a wide range of partial disorder/crystallization 
while still remaining along this characteristic line. These 
runs are the open triangles and squares in Fig. and 
inset. The system can approach arbitrarily close to the 
fully crystallized state before disordering (dilating), but 
can also permanently end up in the fully crystallized 
state. The system displays similar intermittency when 
started from an over-dilated state and allowed to evolve 
to a steady state: during net compaction, the system 
routinely re-dilates slightly before compacting further. 

It is remarkable that at lower 7 the system sustains 
higher forces, in contrast to the classical mechanical or 
thermodynamic dP/dV < 0. This effect is a natural 
consequence of Reynolds dilatancy: shearing causes an 
increase in V. One possibility for the corresponding in- 
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FIG. 6: Dependence of P on F for (a) F = 2.0 and varying 
Q and (b) Q — 0.051 Hz and varying F. Solid symbols are 
crystallized states. Dashed line is fit to (b), shown also on (a) 
for comparison. Inset: Intermittency in run marked A in (b). 

crease in P is that the larger forces are a sign of well- 
developed force chains, unmeasured in this experiment. 
As such, changes in the internal granular pressure would 
depend not only upon the volume, but also upon an order 
parameter, X, describing the strength of the force chains 
present. Thus we could write 

dV ~dV^ \dx)ydV' 

Setting a smaller initial V produces a higher P, indi- 
cating that dP/dV is negative as for ordinary materi- 
als. In 2D sheared systems at constant volume, such as 
shearing increased the strength of the force chains 
and hence the pressure, indicating that the quantity 
{dP / dX)y must also be positive. Therefore, dX/dV^ 
while difficult to measure, must also be positive in or- 
der to obtain dP/dV > as observed in the experiment. 
This formalism suggests two possible interpretations for 
the trend in Fig.Et, where for high Vt dP/dV is positive 
but decreasing. First, there may be an upper limit to 
X due to the number of particles available to participate 
in force chains. Second, if V grows too large, the lateral 
stability of force chains is compromised. 

In this Letter, we have described competing disorder- 
ing and crystallizing effects for sheared/ vibrated granular 
materials. While the concept of a granular temperature 
has provided a thermodynamic analog for granular gases, 
the findings here point to several caveats regarding that 



association. External vibration fails to provide "heat- 
ing" in a thermodynamic sense: the resulting fluctua- 
tions do not allow the shear profile to be liquid-like, and 
they produce crystallization rather than melted disorder. 
However, fluctuations introduced by the shearing of the 
material have the opposite effect. Future work should be 
undertaken to examine temperature-like variables which 
might control this phase transition. 

Granular materials are often studied with constant- 
volume constraints, where they exhibit higher forces for 
higher packing fractions. However, the relationship be- 
tween P and V appears to involve not only the classical 
dP/dV < 0, but a second, opposing, effect dependent 
upon the internal structure. As a result, lower 7 states 
exhibit larger forces. dP/dV is found to have a nearly 
universal value for a wide range of conditions, suggesting 
that predictions for its value should be possible. In addi- 
tion, there is an unexplained logarithmic dependence in 
H(^) for > Qc- 
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